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Abstract. The development of accurate and fast algorithms for the Boltzmann 
collision integral and their analysis represent a challenging problem in scientific 
computing and numerical analysis. Recently, several works were devoted to the 
derivation of spectrally accurate schemes for the Boltzmann equation, but very few 
of them were concerned with the stability analysis of the method. In particular 
there was no result of stability except when the method is modified in order to en- 
force the posivity preservation, which destroys the spectral accuracy. In this paper 
we propose a new method to study the stability of homogeneous Boltzmann equa- 
tions perturbed by smoothed balanced operators which do not preserve positivity 
of the distribution. This method takes advantage of the "spreading" property of 
the collision, together with estimates on regularity and entropy production. As an 
application we prove stability and convergence of spectral methods for the Boltz- 
mann equation, when the discretization parameter is large enough (with explicit 
bound). 
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1. Introduction 

This work is the sequel of several papers devoted to the approximation of the 
Boltzmann equation using fast spectral methods [261 [Iffl EO, US] . The present paper 
is devoted to the stability and convergence analysis of general spectral algorithms. 

In a microscopic description of rarefied gas, the particles move by a constant 
velocity until they undergo binary collisions. In statistical physics, the properties 
of the gas are described by a density function in phase space, f{t,x,v), called the 
distribution function, which gives the fraction of particles per unit volume in phase 
space at time t. The distribution function satisfies the Boltzmann equation, a non- 
linear integro-differential equation, which describes the combined effect of the free 
flow and binary collisions between the particles. 

The main difficulty in the approximation of the Boltzmann equation is due to the 
multidimensional structure of the collisional integral, since the integration runs on a 
highly-dimensional unfiat manifold. In addition the numerical integration requires 
great care since the collision integral is at the basis of the macroscopic properties 
of the equation. Further difficulties are represented by the presence of stiffness, like 
the case of small mean free path or the case of large velocities [9] . 

For such reasons realistic numerical simulations are based on Monte-Carlo tech- 
niques. The most famous examples are the Direct Simulation Monte-Carlo (DSMC) 
methods by Bird [T] and by Nanbu pi]. These methods guarantee efficiency and 
preservation of the main physical properties. However, avoiding statistical fluctua- 
tions in the results becomes extremely expensive in presence of non-stationary flows 
or close to continuum regimes. 

Among deterministic approximations. Discrete Velocity Models (DVM) are based 
on a cartesian grid in velocity and on a discrete collision mechanism on the points 
of the grid that preserves the main physical properties. Unfortunately DVM are not 
competitive with Monte-Carlo methods in terms of computational cost and their 
accuracy seems to be less than first order [221 [231 [211 [31] ■ 

Another class of numerical methods, based on the use of spectral techniques in 
the velocity space, has been developed. The methods were first derived in [25], 
inspired from spectral methods in fluid mechanics [6J and by previous works on 
the use of Fourier transform techniques for the Boltzmann equation [^. They are 
based on approximating in the velocity space the distribution function by a periodic 
function, and on its representation by Fourier series. The spectral method has 
been further developed in f2E[ [28j where evolution equations for the Fourier modes 
were explicitly derived and spectral accuracy of the method was proven. Strictly 
speaking these methods are not conservative, since they only preserve mass, whereas 
momentum and energy are approximated with spectral accuracy. Moreover, the 
spectral method has been applied also to non homogeneous situations [lOj, to the 
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Landau equation [HI EZ] , where fast algorithms can be readily derived, and to the 
case of granular gases [T2]. Independantly A. Bobylev & S. Rjasanow [31 H] have 
also constructed fast algorithms based on a Fourier transform approximation of the 
distribution function, but the method is not spectrally accurate (only second order). 

In [20] a fast spectral method was proposed for a class of particle interactions 
including pseudo-Maxwell molecules in dimension 2 and, most importantly, hard 
spheres in dimension 3, on the basis of the previous spectral method together with 
a suitable semi-discretization of the collision operator. This method permits to 
reduce the computational cost from Oiv?) to 0(n logg n) without loosing the spectral 
accuracy, thus making the method competitive with Monte-Carlo. 

However an important drawback of the spectral methods up to now had been 
the lack of proof of stability. Indeed as compared to discrete velocity methods the 
difficulties are somehow opposite: consistency results are easily obtained, whereas 
the lack of positivity preservation of the scheme is a major issue when one studies 
its stability properties. The only paper concerned with the issue of stability for 
spectral methods applied to the Boltzmann collision operator is [28], but in the 
latter the author introduce some "filters" on the Fourier modes in order to restore 
the posivity- preservation of the scheme, which breaks the spectral accuracy. 

In this paper we give the first stability result for the spectral methods applied to 
the Boltzmann collision operator. Moreover we propose a method which is likely to 
have other utilizations in coUisional kinetic theory: 

• we write the Galerkin approximation on the first N Fourier modes of the 
evolution equation as a a smooth balanced perturbation of the original equa- 
tion, in the sense of a perturbation by some small and mass-preservation 
(although not positivity-preserving) error term; 

• we prove existence and uniqueness of smooth solution for small times, con- 
ditionally to a bound on the norm; 

• we use the mixing structure [29l [T8] of the collision process to show appear- 
ance of positivity after a small time (depending on the size of the box of 
truncation and the approximation parameter N)\ 

• we use the mass conservation to deduce uniform bounds on the L} norm, 
and therefore regularity bounds growing at most exponentially in time; 

• we perform a detailed analysis of the unperturbed truncated problem, show- 
ing uniform in time regularity and asymptotic convergence to equilibrium; 

• finally we use that the equilibrium is unchanged by the smooth balanced 
perturbation, and that it is non-linearly stable for the perturbed periodized 
Boltzmann equation, in order to prove global in time stability and conver- 
gence to equilibrium for the perturbed Boltzmann equation (we connect the 
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previous point for initial times together with the stabihty of equihbrium for 
asymptotic times). 

Hence our paper introduce a general method on how to exploit fine mixing proper- 
ties of the collision process in the study of stability of a particular class of perturbed 
Boltzmann equation, with the application in mind to the stability of spectral meth- 
ods. 

The outline of this paper is as follows. The Boltzmann equation and its basic 
features are presented in Section [2l In Section [3] we explain the truncation and 
periodization associated with spectral methods and fast spectral methods and we 
formulate the problem of stability of these methods in the general framework of the 
stability properties of the Boltzmann equation with respect to a smooth balanced 
perturbation. Section H] is devoted to the proof of the main stability result in the 
general framework. Section [5] is devoted to the study of the asymptotic behavior 
of the truncated problem on the basis of the entropy production theory. Finally 
in Section [6] we apply the latter result to the spectral method and establish some 
stability and convergence results of the numerical solution. 



The Boltzmann equation describes the behavior of a dilute gas of particles when 
the only interactions taken into account are binary elastic collisions. It reads for 
X E ^, V G where 1] G is the spatial domain {d > 2): 



where / := f(t,x,v) is the time-dependent particles distribution function in the 
phase space. The Boltzmann collision operator Q is a quadratic operator local in 
{t,x). The time t and position x act only as parameters in Q and therefore will be 
omitted in its description 



We used the shorthand / = /(f), = /(f*), /' = f{v'), = f{v'^). The velocities 
of the colliding pairs (f , f *) and (f ', v'^) are related by 



2. The Boltzmann equation 



dl 
dt 



+ v-VJ = Qif,f) 



(2.1) 




f^f) dadv^. 
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with a G S"'^^. The colhsion kernel i? is a non- negative function which by physical 
arguments of invariance only depends on 1^ — and cos 6' = u ■ a (where u = 
{v — v^)/\v — Vi,\ is the normalized relative velocity). In this work we are concerned 
with short-range interaction models. More precisely we assume that B is locally 
integrable. Here are the hypothesis on the collision kernel: 

(2.2) B(\u\,cose) ^mu\)b(cose), 
with 

(2.3) $(z) =z^, z e R+, for some 7 G (0, +00) 
and b smooth such that 

/»7r 

(2.4) / 6(cos 9) sin'^-^ 9 d9 < +00. 

These assumptions are satisfied for the so-called hard spheres model B = \v — v^\, 
and it is known as Grad's angular cutoff assumption whea. it is (artificially) extended 
to interactions deriving from a power-law potentials. As an important benchmark 
model for the numerical simulation we therefore introduce the so-called variable hard 
spheres model (VHS), which writes 

B{\u\,cos9) — Cry \u\'^, 

for some 7 e (0, 1] and a constant > 0. 

For this class of models, one can spht the colhsion operator as 

W,/) = Q^iLf) - L{f)f, 

with 

Q^{fJ)= f f B{\v-v,\,cos9)ff:dadv,, 

and 

L(/)=/ / B{\v-v,\,cos9)Udadv,. 

Boltzmann's collision operator has the fundamental properties of conserving mass, 
momentum and energy: at the formal level 

/ QifJ)Hv)dv = Q, cp{v) = l,v,\v\^ 
and it satisfies well-known Boltzmann's H theorem 

/ flogfdv = -[ Q{f,f)\ogif)dv>0. 
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The functional — j f log / is the entropy of the solution. Boltzmann's H theorem 
implies that any equilibrium distribution function, i.e., any function which is a 
maximum of the entropy, has the form of a locally Maxwellian distribution 

/I 12' 

I \u — v\ 

M{p,u,T){v) = —^-rj^exp' ' ' 



(27rT)'^/2 ^^^^ V 2T 

where p, u, T are the density, macroscopic velocity and temperature of the gas, de- 
fined by 

P = [ f{v)dv, u = - f vf{v)dv, T= — f \u — v\'^ f{v) dv. 

For further details on the physical background and derivation of the Boltzmann 
equation we refer to Cercignani, lUner, Pulvirenti p] and Villani [30] . 

3. Formulation of a general stability result 

In this section we remind the basic principles leading to the periodized truncations 
of the Boltzmann collision operator arising in spectral methods. Then, we present 
the main result of this paper: the stability of the spatially homogeneous Boltz- 
mann equation with respect to a smooth balanced perturbation, preserving mass 
and smoothness but not nonnegativity of the solution. This stability means that we 
are able to construct global solutions and estimate the error between perturbed and 
unperturbed solutions. 

Any deterministic numerical method requires to work on a bounded velocity space. 
This therefore supposes a non physical truncation (associated with limit conditions) 
of this velocity space, which we shall discuss below. 

3.1. General framework. We consider the spatially homogeneous Boltzmann equa- 
tion written in the following general form 

(3.1) % = Q{f,f), 

where Q{f, f) is given by 

(3.2) Q{f,f) = I^B{y,z)[ff:-fj]dydz, v eR' 
with 

v' = v + Q'{y,z), v'^ = v + Q'^{y,z), = v + Q^{y, z). 

In the equations above, C is some given (unbounded) domain for y, z, and 9, O', 
0'^ are suitable functions, to be defined later. This general framework emphasizes 
the translation invariance property of the collision operator, which is crucial for the 
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spectral methods. We will be more precise in the next paragraphs for some changes 
of variables allowing to reduce the classical operator fl2.1l) to the form fl3.2p . 

A problem associated with deterministic methods which use a fixed discretization 
in the velocity domain is that the velocity space is approximated by a finite region. 
Physically the domain for the velocity is W^, and the property of having compact 
support is not preserved by the collision operator. In general the collision process 
indeed spreads the support by a factor \/2 in the elastic case (see [29l [18] and also 
[T7] for similar properties in the inelastic case). As a consequence, for the continuous 
equation in time, the function / is immediately positive in the whole domain M*^. 
Thus, at the numerical level, some non physical condition has to be imposed to keep 
the support uniformly bounded. In order to do this there are two main strategies: 

• One can remove the physical binary collisions that will lead outside the 
bounded velocity domains. This means a possible increase of the number of 
local invariants, i.e., the functions (f such that 

(<^1 + - - 

is zero everywhere on the domain. If this is done properly {i.e., without 
removing too many collisions), the scheme remains conservative and without 
spurious invariants. However, this truncation breaks down the convolution- 
like structure of the collision operator, which requires the translation invari- 
ance in velocity. Indeed the modified collision kernel depends on v through 
the boundary conditions. This truncation is the starting point of most 
schemes based on Discrete Velocity Models. 

• One can add some non physical binary collisions by periodizing the function 
and the collision operator. This implies the loss of some local invariants 
(some non physical collisions are added). Thus the scheme is not conserva- 
tive anymore, although it still preserves the mass if the periodization is done 
carefully. However in this way the structural properties of the collision oper- 
ator are maintained and thus they can be exploited to derive fast algorithms. 
This periodization is the basis of spectral methods. 

Therefore, we consider the space homogeneous Boltzmann equation in a bounded 
domain in velocity = [—L,LY {0 < L < oo). We truncate the integration in y 
and z in (13.21) since periodization would yield infinite result if not: we set y and z 
to belong to some truncated domain Cr C C (the parameter R refers to its size and 
will be defined later). For a compactly supported function with support included in 
Bs, the ball centered at with radius S > 0, one has to prescribe suitable relations 
(depending on the precise change of variable and truncation chosen) between S, R 
and L in order to retain all possible collisions and at the same time prevent inter- 
sections of the regions where / is different from zero (this is the so-called dealiasing 
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condition) . Then the truncated colhsion operator reads 



(3.3) Q^(/,/) 



/ B{y,z){rj' - Uf)dydz 
Jcr 



for V G (the expression for f G M'' is deduced by periodization). By making 
some changes of variable on f , one can easily prove for the two choices of variables 
y,z of the next subsections, that for any function ip periodic on Vl the following 
weak form is satisfied: 



(3.4) 



/ Q^{fJ)^{v)dv = ^f [ B{y,z) f^f {^'^ + <p' -<p^-^) dydzdv. 
Jvl ^ Jvl Jcr 



Now, we use the representation to derive spectral methods. Hereafter, we 
use just one index to denote the ci-dimensional sums with respect to the vector 
k = {ki, .., kd) E Z*^, hence we set 

N N 

E ■ 

k=-N ki,...,ka=-N 

The approximate function is represented as the truncated Fourier series 



N 



(3.5) Mv) = J] A 

k=-N 

with 



{2LY Jj,/^ ' 

In a Fourier- Galerkin method the fundamental unknowns are the coefficients /^(t), 
k = —N, . . . ,N. We obtain a set of ODEs for the coefficients fk by requiring that 
the residual of (13. 3p be orthogonal to all trigonometric polynomials of degree less 
than A^. Hence for k = —N, . . . ,N 



(3.6) ^^-QR(f^J^)je-^-Lk-^dv = 0. 

By substituting expression (13. 5p in (13.41) we get 

Q^'ifN, In) = Q^'+(/^, In) - L'^iM h 



ANALYSIS OF SPECTRAL METHODS FOR THE BOLTZMANN EQUATION 



with 

(3.7) L^{fN)fN = E E Pirn,m)fJ^e^i^'^"^>\ 



N N 



l=-N m=-N 
N N 

{l+m)-v 



(3.8) Q^'+(/^,/^) = E E Pil^rn)fJme'T 

l=-N m=-N 

where 

(3.9) p(l,m)= [ B{y,z)e'T^i'-®'^y''^+"'-®'*^y''^) dydz. 



JCr 

The spectral equation is the projection of the colhsion equation in P^r, the (2A^ + 
1) '^-dimensional vector space of trigonometric polynomials of degree at most in 
each direction, i.e., 

= Q^Un, In), 
where Vn denotes the orthogonal projection on P^r in L'^(T>l). 

3.2. The truncation associated with classical spectral methods. In the clas- 
sical spectral method a simple change of variables in fl2.ll) permits to write 



(3.10) Q(/,/)= / 


/ B^^^^^{g,uj){f{v')f{v: 






with g = V — E M'^, to E ! 


3*^^^ and 




( v' = v- \{g - \g\uj), 


(3.11) 


I v'^ = v -\{g+\g\uj), 




y v^ = v + g. 



Then, we set C := M"^ x and 

^'{g,^) ■= —{g ~\9\^), Q'^ig,^) ■= -^{g + \g\^), Q*ig,^)-=g- 

Finally the collision kernel B'^''"''^^ is defined by 

(3.12) B^'^'^ig, CO) = 2""' {l - (g ■ co)Y^'~' B{\g\,2{g ■ cof - l). 

Thus, the Boltzmann operator (13.101) is now written in the form (13.21) . We consider 
the bounded domain Vl = [—L, L^, (0 < L < oo) for the distribution /, and the 
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bounded domain Cr = Bji x S'^"^ for some R > 0. The truncated operator reads in 
this case 

(3.13) Q^{fJ){v)= [ B^^'^^^{g,u){f{v:)f{v')-f{v.)f{v))dudg. 

3.3. The truncation associated with fast spectral methods. Here we shall 
approximate the collision operator starting from a representation which conserves 
more symmetries of the collision operator when one truncates it in a bounded do- 
main. This representation was used in [3l [T^ to derive finite differences schemes 
and it is close to the classical Carleman representation (cf. [7]). The basic identity 
we shall need is (for u G Mf^) 

(3.14) I [ F{\u\a -u)da = [ 5{2yu+ \y\') F{y) dy. 
Using (13.141) the collision operator (12.11) can be written as 



(3.15) Qif, f){v) = 2"'' / / B^'^-'iy, z) 5{y ■ z) 

{f{v + z)f{v + y)- f{v + y + z)f{v)) dy dz, 

with 

B^'^^'iy, z) = 2'^' B(\y + z\, -^-^^±41 \y + z\-^''-'\ 

\ \y\\y + z\J 

Thus, the collision operator is now written in the form (13. 2p with C := M'^ x W^, 

Biy,z)=Bf^'\y,z)6iyz), 

and 

v[ = V + Q'^{y,z), v' = V + Q'{y,z), = v + <d^{y, z). 

with 

0Uy,^):=^, Q'{y,z):=y, e^{y, z) := y + z. 
Now we consider the bounded domain = [—L, L^, (0 < L < oo) for the 
distribution /, and the bounded domain Cr = Br x Br for some i? > 0. The 
(truncated) operator now reads 
(3.16) 

Q\f. f) (v) = I B^^^^'iy, z) 5(y ■ z) {f{v + z)f{v + y)-f{v + y + z)f{v)) dy dz, 

JBrXBr 

for V G Pl- This representation of the collision kernel yields better decoupling 
properties between the arguments of the operator and allows to lower significantly 
the computation cost of the method by using the fast Fourier transform (see [20} [13]). 
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Let us make a crucial remark about the choice of R. When / has support included 
in Bs, S' > 0, it is usual (see [221 ED]) to search for the minimal period L (in order 
to minimize the computational cost) which prevents interactions between different 
periods of / during one collision process. From now on, we shall always assume that 
we can take L and R large enough such that, when needed, R > \/2 L. Hence all 
the torus is covered (at least once) in the integration of the variables {g, uj) or {y, z). 

3.4. A common abstract formulation for the stability of spectral methods. 

From now on, shall denote a periodized truncated collision operator as in f l3.13p 
or (13.161) . As we shall see, using this formulations, both classical and fast spectral 
methods fall into the following framework: 

{ % = Q''{fJ) + Peif), veVr., 
(3.17) <^ 

where is a "smooth balanced perturbation", which means that it satisfies the 
following (balanced law) 



(3.18) / Peif)dv = 

and preserves the smothness of the distribution function, i.e., there exists constants 
Co, Cfe > such that 

r \\Pe{f)\w < Coii/iu^ii/iui 

(3.19) 

i lin(/)lk,.„, < fc>0, 

where || ■ H^fc^^ is the usual norm of the Sobolev space of periodic functions H^^^iVi). 

Moreover the perturbation is supposed to be small in the following sense: there 
exists a function ^{e) such that for any p > 0, 

(3.20) \\Pe{f)\\Hl.^ < vie), 

where ip{e) depends on ||/||^p+fc for some A; > 0, and goes to zero as e goes to zero. 

Finally in order to prove global existence with uniform regularity bounds, we 
shall require additionnal assumptions on the relation between the equilibrium dis- 
tributions of the perturbed and unperturbed (periodized) Boltzmann equations, and 
about the stability of the unperturbed equation (see the following statement). 

Let us therefore write the unperturbed equation for reference: 

(3.21) { l = '3"(^-^)- 
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Let US state the general stability theorem: 



Theorem 3.1. Let us consider a perturbed Boltzmann equation i3.11\) in the torus 



T?L, where is defined by Ii3.13\) or Ii3.16\) . and for a sequence of smooth balanced 
perturbations {P^ = Peif))^^^ which satisfy (l338D - (l339D - fl3:20D . 

Assume that the constant functions are equilibra of the pertubed equation fl3.17p 
(as for equation (13.211) ) and that they are nonlinearly locally stable in any H^^^{T>l) 
for equation f l3.2ip . 

We assume that /o is a non-negative function, non zero everywhere, belonging to 
Hperi'^L) with A; G N and k > d/2. We consider a sequence of smooth balanced per- 
turbations /o,£ of the initial datum for the perturbed problem ( 13.171) (non necessarily 
positive) such that 

fo,e= / /o and ||/o - /c^ll^fc < 

with ip{e) goes to zero when e goes to zero. 

Then, there exists eo > depending only on the collision kernel B, the truncation 
R, the constants in fl3.19l) - fl3.20p for the perturbation, and the L^(T>l) and H^^^{T>l) 
norms on fo, such that for any e G (0, ^o); 

(z) there exists a unique global smooth solution f^ to fl3.17p ." 
{ii) for any p < k, this solution belongs to H]p^^{T>i) for all times with uniform 
bounds as time goes to infinity; 

(iii) this solution remains "essentially non-negative" uniformly in time, in the 
sense that there is rj{e) > (with r]{e) as e goes to 0) such that the 
non-positive part is r]{e) -small: 

vt>o, \\fr{t,-)\\L^<v{s) 

where f~ denotes \fe\ l{/,<o}; 

(iv) this solution f^ converges in H^^j.{T>l) for any p < k, uniformly on any 
[0,T], T > 0, to the solution f of the unperturbed equation f l3.2ip when the 
parameter e goes to zero; 

{y) the solution fe to fl3.17p converges in Hp^^{T>L) as time goes to infinity to 
the constant equilibrium distribution in the torus prescribed by its mass, and 
it is "asymptotically uniformly positive" , that is for t larger than some fixed 
explicit time. 

We split the proof into two main steps: first in Section H] we prove existence, 
uniqueness and smoothness of a solution on an arbitrary bounded time interval (as 
the size of perturbation goes to 0). The main difficulty is to prove that non- negativity 
of the distribution function is recovered in a certain sense. Then, in Section [5] we 
study the asymtotic behavior and establish global stability in time. The main issue 
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is there to prove regularizing properties of the gain operator Q^'^{f, /) and entropy 
production estimates on Q^. Finally in Section El we apply the previous general 
results to spectral methods, and prove their stability and convergence. 

4. Proof of stability on an arbitrary bounded time interval 

In this section we first give some technical lemmas and next establish a result 
showing existence and uniqueness of a smooth solution on an arbitrary time interval 
to the perturbed equation fl3.17p . depending on an assumption of smallness on the 
size of the perturbation. Then, in Lemma 1^31 we prove the control of negative values 
of fit). 

4.1. Preliminary results. We start this section by a classical result of estimates 
on the Boltzmann operator Q^{g, h) given by 

Q\g,h):= [ B{y,z){g:h' - g,h)dydz. 

JCr 

Lemma 4.1. Let the collision kernel B satisfy the assumption f l2.2p -f l273l) - fl2.4p . 
Then, the periodized Boltzmann operator (defined by i\3.13\) or Ii3.16\) ) satisfies: 
for all p G [0, oo] there exists a constant Cp{R, B) > such that 

(4.1) \\Q''ig,h)\\L., WQ^'iKg^L. < Cp{R,B)\\g\\L^\\h\\L,. 

Proof. The proof is exactly similar to the case of the usual Boltzmann collision 
opertor for a collision kernel bounded with compact support, see for instance [1F| 
for a recent proof. □ 

Now, we prove smoothness of the solution to the perturbed problem ( 13.170 on a 
fixed time interval under the assumption of an a priori bound on the norm of 
the solution. 

Lemma 4.2. Let us consider a collision kernel B which satisfies the assumptions 
fl2.2p - (12.30 - 02. 4p and a sequence of smooth balanced perturbations (P^ = Pe{f'^))^yQ 
which satisfy (13. 180 - (13. 190 - (13.200 . and let T > be the length of the time interval. 
Assume that /o G H^{Vl) for k en and that f{t) is a (non necessarily positive) 
solution to (I3.17P with initial datum fo, which satisfies the L^-estimate 

(4.2) vtG[o,r], ||/(t)|Ui<M. 

Then, there exists a constant Ck{M) > 0, only depending on M, R, T and ||/o||_fffc^,, 
such that 



(4.3) 



vtG[o,r], Wfmn^^^ < Ck{M). 



14 FRANCIS FILBET AND CLEMENT MOUHOT 



Proof of Lemma 4-2 We proceed by induction on k >0. For the first stage k = 0, 
we apply Lemma 14.11 with p = 2 and g = h = f and we use assumption (13.191) on 
the perturbation: 

< {C,{R,B)+Co) ||/(t)|U.||/(t)||i.. 

Hypothesis (14. 2p provides a control on ||/(t)||Li and we can apply Gronwall's lemma 
to get the result (14.31) at stage k = 0. 

Let us now assume that (14. 3 p holds at stage k > and let us prove that it also 
holds at stage (/c + 1). 

Let us first recall a formula on the derivatives of Q^: from the bilinearity of 
and the translation invariance property of the periodized Boltzmann collision 
operator, one has 

V.Q^(/, /) = Q^(V./, /) + Q^(/, V./), 
which yields a Leibniz formula at any order s G N: 



\u\<s 



(4-4) < E (l)m9^f^9^-'mi- 

Now, using (14. 4p with s = k + 1 we have 
(4.5) 

mfj)\\i^.^ < cmfj)\\%^^ + c y: e ( ' ) md^f^d^-^m 



2 

L2- 



\u\=k+l |Ml<fc+l 



From Lemma 14.11 with p = 2 and g = d'^f, h = d'^ ^ f together with the hypothe- 
sis (14. 3p . we get 

(4.6) ||Q^(/,/)L.^^ < C(A:,i?)||/||^.^J|/||^.^^ < C,{k,R,B)C,{Mr. 

Then we split the last term of (14.50 in two parts for /i 7^ and fi = 0. For yU 7^ 0, we 
again apply Lemma [4.11 with p = 2 and g = d^f, h = d'^~^f and use the fact that 
both derivatives = XliLi 1/^*1 — ^ ^^"^ \^ ~ — ^- 

(4.7) J2 E ( u ) \\Q''id''f,d-'^f)\\l. < C,{k,R,B)C,{Mf. 

\v\=k+l lMl<fc+i ^ ^ 
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Finally, for fi = we apply Lemma [4.11 with p = 2 and g = f, h = d'^f: 

(4.8) <C,{k,R,B)\\f\\,.\\f\\^.,.. 

\u\=k+l 

Then, gathering inequalities from (14. 6 p to (14.81) and using the assumption (I3.19P on 
the smooth balanced perturbation P^, we have 

< \mfj) + Ps{f)h^.,\\fh^.^ 

< C{C2 + C,) Cu{Mf 11/11 ... + C (C4 + C) ll/IUi 11/115,..!. 

-•-'per iip^j- 

Finally using the control (14. 2 p on ||/(t)||Li we apply Gronwall's lemma to get (14. 3 p 
at stage k + 1: there exists a constant Ck+i{M), only depending on M, R, T and 
ll/oll<+i such that 

VtG[0,T], ||/(t)||^.+i < C,+i(M). 

□ 

Then, we establish existence and uniqueness of a smooth solution for the perturbed 
problem (I3.17P on a small time interval [0, f], f > 0. 

Proposition 4.3. Let us consider a collision kernel B which satisfies the assump- 
tions (I2.2p - (l2.3p - (l2.4p and a sequence of smooth balanced perturbations (P^ = PeiD^yQ 
which satisfy fl3A8|) - (l3A9D - (13:201) . We assume that fo G H^^^{Vl), forkeN and 
set 

(4.9) M = 2||/o|Ui. 

Then, there exists f = f{M) > such that for all e > the perturbed Boltzmann 
equation (13.17^ admits a unique solution ( non necessarily positive ) on the time in- 
terval [0,f], where the solution f satisfies 

(4.10) VtG [0,f], ||/(t)|Ui < M. 

Moreover, there exists a constant Ck{M) > 0, only depending on M, R, and ||/o||_fffc^^ 
such that 

(4.11) VtG[0,f], ||/(t)||^.^^ < Cfc(M). 

Proof of Proposition First, we apply Lemma [4. II with p = 1 and g = h = f 

||qV,/)||,, <Ci(P,p) ll/iu. ll/iu.. 

Moreover, using assumption (I3.19p on the perturbation, there exists C > such that 
for all £ > 

\\Ps{f)\\L^<C\\f\\h. 
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Therefore, we obtain a constant C > 0, only depending on R and the colhsion kernel 
B such that 

j^\\m\\L^<c\\f\\i. 

This implies that 

J- - L II/oIIli * 

Now, setting M = 2 ||/o||li and from the latter inequality, we show that there exists 
f < l/(2C||/o||ii) such that 

Vt e [0,f], ||/(t)|Ui<M, 

which gives (14. 101) . 

From the estimate in L}{T>l) on the function /(t) on the time interval [0,f], 
we prove existence and uniqueness of a solution by Cauchy-Lipschitz theorem in 
L}(T>l) (because of the truncation on the collision kernel is a bounded bilinear 
function from L^{Vl) x L^{Vl) to L^{Vl)). Finally, from the bound KW^ and the 
smoothness assumption /o G H^^^.{Vi) on the initial datum, we are able to apply 
Lemma [4.2^ which proves that there exists a constant Ck{M) > 0, only depending 
on M, R, T and ||/o||_H-fc such that 



Wfrnn^^^ < Cu{M). 

This concludes the proof. □ 

By iterating Proposition 14.31 and Lemma 14. 2[ we observe that uniform control on 
the L}{T>l) norm on an arbitrarily large time interval [0, T] together with smoothness 
on the initial datum will ensure existence and uniqueness of a smooth solution on 
this time interval [0,T]. Furthermore we observe that the control on the L}{T>i) 
norm is obvious for the classical Boltzmann equation thanks to the positivity and 
mass preservations. Therefore, we shall now focus on the control of positivity of the 
solution, showing that the solution remains "almost positive" for arbitrarily large 
time interval, hence allowing to produce uniform control on the L^{T>l) norm for 
arbitrarily large time interval. 

We first state a technical lemma which takes advantage of the mixing property of 
[yi order to show spreading of the support of a characteristic function of a ball. 

Lemma 4.4. Let us consider a collision kernel B which satisfies the assumptions 
0-(I23D-(I23D and a truncated operator defined by BAR) or BAR) . Then for 
all < r < V2L, we have 

for some explicit fi = fi{R, L) > 1 and Cq. 
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Remark: Note that for r > \/2L^ one has ^-B{y,r) = 1 on the torus [—L^LY and 
there is nothing to prove. 



Proof of Lemma 4-4 The invariance by translation allows to reduce the proof to the 
case V = 0. The invariance by rotations implies that Ir := Q^'^{lB{Q,r), ls(o,r)) is 
radially symmetric. More precisely, taking a C°° radially symmetric function (j) such 
that (f) > on B = B{0, r) and < 1b on M.'^, we have 

• the function v — > Q^'^{(f), (f)){v) is continuous, 

• for all V G Ir{v) > g+(0,0)(f), 

• for all veB, g^'+(0, 4i){v) > 0. 

As a consequence, for any ball B' = B{0,r') strictly included in B, there exists 

Kr' > such that Ir > Kr' 1 B' ■ 

In order to conclude, we just need to estimate the support of Ir close to the ball 

B. 



Let us fix r' G (0,r) and choose v',v'_^ G -8(0, r') such that \v' 
min{i?; \/2r'}. Then taking a in -8(0, r') such that 
we find \v\ = \/2r' if ^/2r' < R and 



\v'-v[ 



r' and 



f , f^, f', is a square. 



p2 

[r 



1/2 



R 
+ 2 



else. We define 

/io = min 1 — 

R/V2<r'<V2L 



R^ 



4(r')2_ 



1/2 



R 

H = mm 

2r H/(2^/2L)<1/<l/^/2 



l-y^ + y]e (1,V2). 



This concludes the proof: we deduce that for any v such that l^;! G (r',/ior'), we 
have 

Ir{v)>Q+{<PA){v)>^ 

since is strictly positive in the neighborhood of the f ', v'^ associated to v constructed 
above. Hence we deduce by taking r' < r close to r that for any < fi < fiQ, we 
have 

Ir > C{fl,r) lB(0,/ir) 

for some constant C{fi,r) > depending continuously on r. We can choose /i = 
(1 + /xo)/2 for instance, and, for this choice of /i, we take 



Co = min C(yU, r) > 0. 

0<r<V2 



□ 



Finally we establish the following positivity result on the solution to the perturbed 
problem fl3.17p : 



18 



FRANCIS FILBET AND CLEMENT MOUHOT 



Lemma 4.5. Let us consider a collision kernel B which satisfies the assumptions 
(12. 2p -( 12731) -f l2.4p . a truncated operator defined by Ii3.13\) or Ii3.16\) . and a sequence 
of smooth balanced perturbations (P^ = Pe{f))^^Q which satisfy (13 . 1 8p - ( l3?T9l) - (13 . 201) . 

We assume that fo is a non-negative function such that fo G H'^^^{'Dl) with k G 
N and k > d/2. Moreover, we define M = 2||/o||li and fo^s o smooth balanced 
perturbation of fo, which is non necessarily positive and such that 



fo,e= / fo and \\fo - fo,e\\m,r - ^(^) 

where ipi^) 9oes to zero when e goes to zero. We also set f, Ck{M) > the constant 
defined in Lemma \4-S\ such that from Proposition \4-3\ we have 



yte[0,f], Umn^^^ < Cu{M). 



Then, there exists f G (0, r) which only depends on M , R and the collision kernel 
B, and there exists e > which only depends on f, Ck{M), such that for all e such 
that < e < e and for any smooth solution H^^^{T>l) to the perturbed Boltzmann 
equation (I3.17P with perturbed initial datum /o,e, we have 

^v^Vl, f,{f,v)>0. 

Moreover there exists ri{e), which goes to as e goes to zero, such that the non- 
positive part of f satisfies 

(4.12) lirWIU- < vie), tG [0,f], 

Proof of Lemma \4-5\ Let f > be the length of the time interval for which there 
exists a smooth solution to the perturbed Boltzmann equation (13.170 with perturbed 
initial datum /o.e such that (in the following we omit the subscript e for the solution) 

||/(i)||^.^,. <C,(M), t G [0,f]. 

We split the proof into three steps: first, we give a classical estimate on the loss 
term L^{f), second we establish an estimate of /~ with respect to the amplitude of 
the perturbation Pdf), and third we use the spreading properties of Q^'~^{f, f) to 
prove that there exists f G (0, f) such that 

/(f,i;)>0, vEVl. 



Step 1. Applying Proposition 14.31 for the control of ||/(t)||Li on the time interval 

[0,f], we get 

(4.13) \\L\f)\\L^ < C{R,B)\\f\\L^ < C{R,B)M, 
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which gives for all t G [0,ro], with tq = rain {f,\n 2 / (M C{R, B))} 
(4.14) 2 > e*^'^(^-^)* > e--^o^''(/W)^'^ > (,-MC{R,B)t > i_ 



2 



Step 2. Let us split / as / = /+ — / , with = max{0, ±/} and use the 
monotonicity of Q^'^ for nonnegative distribution functions; it follows that 

Q'^'Hfj) = Q'''^{r-r,r-r) 
(4.15) > - [Q''^^{f\f-)+Q''^^f-,n]. 

On the time interval [0,ro], we apply Proposition 14.31 to estimate ||/(t)||Li and since 
< II/WIIli; we get from LemmaO 



(4.16) ||g^'+ (/+,/-) WQ^^'^ir^nh- < C^{R,B)M\\f-{t)U^, 
which yields using (14.151) 

(4.17) Q'^'+l/,/) > -2CMB)M\\rit)\\L^. 

Thus from the Duhamel representation of the solution /, we have for v G V^, 

+ r [g^'+(/(s),/(s)) + P.(/(s))] {v)e~fsL-if(umv)du^^^ 
Jo 

> -2^(e)- r[Q^'+(/(.),/(s)) + P,(/(3))](i;)e-/^'^"(^("»('')'^"cis. 

^0 

Hence, we get from the lower estimate (I4.17P of Q^''^{f,f) and the smallness as- 
sumption (13.201) of the perturbation Pe{f), for all v G Vl 

r{t,v) = max{0, -f{t,v)} < 2^{6)+2 [ {2CUR, B) M \\f- {s)\\l^ + ^{e)) ds. 

Jo 

Finally, we take the supremum in f G and apply Gronwall's lemma to get for 

any < r < Tq 

lirWIU^ < (2^(e) + 2ro^(e))e^'^-(^''^)*^*:=r^(5), te[0,T,], 

which proves (14.121) on the time interval [0,ro]. 

Step 3. Let us prove that there exists f G (0, tq) such that 

f{f,v)>0, vEVl. 
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We start again with the Duhamel representation of the solution 



but we now take into account the fact that the first term is essentially positive and 
we use the spreading property of the operator Q^'^ (Lemma I4.4p . 

On the one hand since the initial datum is smooth enough {k > d/2 is large 
enough such that /o,e is Holder), there exists an explicit 6 > depending on Cfc(M) 
such that for e small enough, there exists Vq eV^ such that 

/o,e(f) > |lij(t,o,5)(t^) - V'(^), withr; = 1^^- 

On the other hand, using the lower bound on the gain operator fl4.17p and the 
estimate (14.121) of /~(t), it gives 

Q'^'^ifJ) > -2CUR,B)Mr]{e) 
Finally, using (14.141) and the smallness assumption on Ps{f), it first yields for any 

< T < To 

f{t,v) > AolB{vo,5)iv) - (pi{e), 

with 

(4.18) A:=|, M^)-='^ro{2C^{R,B)Mr]{e) + ^ie)) + ^{e). 

Now we choose £i > small enough such that 

< ^ < Ao - v^o(^i) ^ Mei) < Y- 



Thus, we get on any time intervall [0, r] C [0, tq] 



Hence, using the spreading properties of the operator Q^'~^{f, f) of Lemma [4.41 and 
the monotonicity of Q^'~^ for nonnegative distribution functions, it follows that 



42 

4 



> —r-ColB{vo,fi5)- 
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Next, we again use the uniform bounds previously established in (14 .161) on Q^'^{f^, f~) 
and Q^'+(/-,/+): 

Q^'+(/,/) > g^'+(r,/+)-Q^'+(r,r)-g^'+(r,r) 

> -fColB(voM)-'2C^iR,B)M{ri{e) + M£))- 

Finally, from the smallness assumption fl3.20l) of the perturbation Pe{f), it yields 
for t e[T/2,r] and with fITOD 

f{t,v) > -ij{e)+ f[Q'''^{f{s),f{s)) + Psifis))]{v)e-fsL''^f(-m^^<^-ds, 
Jo 

> ^^ColBivo,,s)-2T [2CMB)M{r]{e) + Me)) + - ^P{e) 

with 

A, = ^AlCo, Me):=2ro [2C^{R, B) M {7]{e) + M^)) +^{e)]-^{e). 

Now, we proceed by induction: assume that there exists {Aj,ej,ipj) such that on 
the time interval [r — t/2^ , r] C [0, Tq] and for e G (0, Ej), we have 

f{t,v) > Aj lB(vo,^^^ 5){v) - ^ji^), 

where 

and v^j(e) — i> as e — 0. 

Using the same method as before, we first set e^+i such that 

and prove that for t E [r — i^+r, t] 

Q'^'+a,/) > (y) ColBM+^8)-2CUR.B)M{r^{e) + ^,{e)). 

Then, from the Duhamel formula and the smallness assumption of Pe{f), we finally 
get the following lower bound 

f{t,v) > Aj+ilB{vo,fi^+^5) - Vj+li^), 

with 

A,+i = JaJCo, ip,+i{e):=2ro [2C^{R, B) M {7^{e) + ^.{e)) +^i6)]-^{e). 
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Since /i > 1, the ball B{vo,^^6) eventually recovers the periodic box [—L,LY 
i.e., for some J large enough: [—L,LY C B{vo, fi"^ S), and for all t E [r — jj,t], by 
applying J's times the previous induction we get for e G (0,ej): 

f{t,v) > AjlB(^vo,,i-'S)iv) - Vji.^)- 

Finally, up to reducing e further, we have proved that there exists (f, e) which 
only depend on the collision kernel B, the initial datum /q, L and the perturbation 
function ip = ip{e) such that for all < e < e, 

\/veVL, fif,v)>0. 

□ 

4.2. Existence and regularity on a bounded time interval. 

Proposition 4.6. Let us consider a fixed time T > 0, a collision kernel B which 
satisfies the assumptions (12.21) -( I2l3l) -f l2.4p . a truncation defined by 1^3. 13\) or 
113.16]) . and a sequence of smooth balanced perturbations (Pg = Pe{f))^~^Q which 
satisfy (IXT8l) - (l3T9l) - (l3:20l) . 

We assume that /o is a non-negative function, not zero everywhere, and such that 
fo e H!^^^{Vl) with k e N and k > d/2. We define M = 2 ||/o||li and {fo,e)e>o a 
sequence of smooth perturbations of fo (which is non necessarily positive) such that 

/ fo,e= / fo and ||/o - /csIIh* < i'i^) 

where ip{e) goes to zero when e goes to zero. 

Then, there exists e > 0, which only depends on the H^^^iVL) and L^{Vl) norms 
of fo, such that for all e G (0, e), 

(i) there is a unique smooth solution f^ = ■) on [0,T] to the perturbed 

equation (13.171) with initial datum /o,^; 
{ii) this belongs to Hp^^iVi) (with bound growing at most exponentially); 
[Hi) there is some explicit r]{e) > (with rj{e) — >• as £ goes to 0) such that the 
non-positive part is ri{e)-small: 

vtG [o,r], ||/-(t,-)IU^ <^(^) 

where f~ denotes max{0, —fe}; 
(iv) this solution satisfies for any p < k 

VtG [0,T], < ^(e) 

where f{t, ■) is the solution of unperturbed periodized Boltzmann equation 
Ii3.21\) . and ip{e) is another explicit function which goes to zero as e goes 
to zero. Hence up to reducing e, the perturbed solution remains close to the 
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unperturbed solution on the finite time interval on which we have constructed 
it. 

Proof. We set 

(4.19) M,{T):=\\U\U.^^e''^^. 

First, applying Proposition 14.31 we have proven that there exists a small f > such 
that the perturbed Boltzmann equation fl3.17p admits an unique smooth solution on 
the time interval [0, f ] with 

mt)\w < M 

and 

Moreover from Lemma 14. 5^ there exist f < f and e > 0, only depending on M, 
Mk{f) < Mk{T), R and the collision kernel B such that for all < e < e, 

\fveVL, fe{f,v)>0, 

lirWIU- < tG [0,f]. 

and 

Then, from the preservation of mass under the action of and P^: 

Q''{feJe){v)dv = [ Pe{fe){v)dv = 0, 

we have that 

/ fe{'r,v)dv = / fo{v)dv. 
Since /o is a non negative function, it gives that at time f 

\\fe{f)\\L^ = ||/o||li, 

and on the time interval t E [0, f] we have 

ll/.WIUi < M and ||/e(t)||^,.^^ < M,(T). 

Therefore, we consider the perturbed Boltzmann equation (13.171) starting from /^(r) 
as initial data. On the time interval [f, 2f], we apply Proposition 14.31 and get that 

||/.(t)||Li < 2||/,(f)|Ui = 2||/o|Ui = M, Vt G [f,2f] 

and 

< mUe'''^ = M,{T), Vt G [f,2f]. 



24 FRANCIS FILBET AND CLEMENT MOUHOT 

Moreover, since f only depends on M, Mk{T), B and i?, we can again apply 
Lemma [4.51 on the time intervall [7^,27"], which yields that 

/.(2f) >0. 

We finally proceed by induction to prove existence and uniqueness of a smooth 
solution /g of the perturbed Boltzmann equation (13.1 7p on the time interval [0,T], 
which proves assertions (i), {ii) and {iii). 

To prove (iv), we compute the difference between the solution f{t) to the unper- 
turbed problem (13.211) and the solution f^{t) to (13.171) : 

^^^^^ = \ {Q\f -feJ + fe) + Q\f + fej- fe)) + Peife)- 

Then, using the smoothness of / and f^, we have from Lemma [4.21 for any p < k 

\\Q\f-feJ + fe)\\H^,,^, WQ^'if + feJ-fenH-,^^ < C,(M) ||/ + || J| / - 

and since the perturbation is small (assumption (I3.20p ) 
it yields that for all t G [0, T] 

for some function (^(e) going to zero as e goes to zero. □ 



5. Asymptotic behavior and global in time stability 

In this section we shall study the asymptotic behavior of the (unperturbed) pe- 
riodized Boltzmann equation (I3.2ip based on a regularity study (in the spirit of 
[T9] ) and the entropy - entropy production theory (mainly relying on the method 
developed in [31]). Finally on the basis of these results we shall prove a global in 
times stability result for the perturbed equation (13.170 . 

5.1. Regularity study of the periodized Boltzmann equation. Let us prove 
the following result 

Proposition 5.1. Let us consider < /o G L^(Vl) such that /q G Hp^^iVi) for 
some k > 0. Then there is a constant C > depending on the and H^^^iVi) 
norms of fo such that the unique global non-negative solution {f{t))t>o to the peri- 
odized equation ^3.21\) satisfies 



Vt>0, ||/(t)||^.^^<C. 
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We proceed as in [12]. In particular we shall extend Lions regularity result on 
to the truncated case Q^'^ [13 HB] . Hence we shall first prove the regularity 
property on the gain operator when the collision kernel is smooth and compactly 
supported, avoiding cancellations at zero relative velocities. Then, we shall include 
the non-smooth part of the kernel using the loss operator. 

We shall split the collision kernel into a smooth and a non-smooth part. As a 
convention, we shall use subscripts "s" for smooth and "ns" for the non-smooth 
parts. In terms of the classical truncation (13.131) we set 

Bf{\g\,cose) = B{\glcose) x^ilgl) U\9\) e,(cos^), 
B^M,oos9) = B^{\g\,cos9) - B^{\g\,cos9), 

where is the mollified version of l|g|<_R, Cnid) is the mollified C°° version 

of l|g|>,;, and is a function on [—1, 1] which is 1 on —1 + 2?7 < m < 1 — 2?7, 
and in [—1, —1 +77) and (1 — 77, 1] (the parameter r} is the mollification parameter). 
In terms of the "fast" truncation (13.161) we set 

Bf{\y\, \z\) = B{\y\, \z\) x^(k|) x^M) C.(k + (^^p 

BU\y\M) = B\\yl\z\) - Bf{\yl\z\l 

with the same notations. 

We deduce the following decomposition of the collision operator: 

where (for instance, with the variables from the "fast" truncation) 
(5.1) Qf'+(/,/)= I [ Bf{\yl\z\)6{yz)f{v + z)f{v + y)dydz. 



Under the assumption that both $ and b defined in (12. 2p are smooth, the regularized 
truncature introduced above ensures that there exist two functions $^ and such 
that 



(5.2) 



B^{\z\, cos 6) = <!>^{\z\)b^{cos 6) 

In the following lemma we shall prove the regularity property of Qf'"*". 

Lemma 5.2. Let B^{\v — t'^|,cos6') satisfy the assumption l{2.^) - (K4\ ) and h5.^) . 

Then, for all r G M"^ 

||Qf'+(/,/)|| .+ -1 < C,e,(r,i3f 



,r+^ ^ ^regVr.O^ ) \\J 
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where the constant Creg{r, Bf) only depends on r and on the collision kernel. 

Proof. We closely follow the proof given by Lions fi5[ [16] and simplified and then 
reformulated in [32l |T9]. Again the preservation of the translation invariance by the 
truncation is fundamental. Starting from the collision operator in the form (15 .ip 
and performing a change of variable we get for v eVl 

Qf '+(/, /) = / BfiK - v'\, \v' -v\)6 (K - v) . iv' - v)) /«) /K) dv: dv\ 

where Bf only depends on Bf. Then we set [33] 

Tg{y)= I Bf{\zl\y\)g{z)dz, T,g{-) = g{- - z), 
Jy+y^ 

where 

y^ = {zeW', z'-y = Q] 

and easily get for t> G 

Qf '+(/,/) = / f{v'){r,,oTor^,,)f(v)dv'. 

Now, we want to estimate Sobolev norms of <5f as a function defined in the torus 
T>L. Applying the Fubini theorem with the discrete and Lebesgue measures and the 
Cauchy-Schwarz inequality, it leads to estimate the Sobolev norms of the Radon 
transform T on the torus T>l 

(5.3) ||gf'+(/,/)||' ^ < ll/IUi / \f{v')\\\T,>oToT^,,f\\\^^dv'. 

On the one hand, since the kernel Bf is compactly supported in y and the 
operator T maps periodic functions g to a. compactly supported function Tg with 
supp {Tg) C Br C Vl. Then, we can consider Tg as a function in the whole space 
. 

On the other hand, using the regularized truncations ^"^^ ©r? and the smooth- 
ness of $ and h in (12^ . it yields that = ^^{\v - v^\) b^{cos9) with 

Then, we can directly apply the result in [l_9j , where the authors proved the following 
regularity estimates on the Radon transform T for smooth kernels 

||T^7||^.,^ < Cre,{r,Bf)\\g\\Hr 

for a function g defined in W^. However in the proof of the latter inequality, we 
can replace g by the smooth and compactly supported function g for which 
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supp {x^ g) C -8(0, R) C Vl. Thus, for all g defined in the torus we get 
(5.4) \\Tg\\,-,i^ = WTgW^.,^ < Creg{r,Bf)\\g\\Hr^^. 

Finally gathering (15.31) and (15. 4p . we obtain the result 

||gf'+(/,/)|| .+-1 < ae,(r,sf)ii/iui||/b,v. 

^per 

□ 

Corollary 5.3. Let B^{\v — v^\, cos 9) satisfy the assumption /12.^) - (K4\ ) and ( 15. 

Then, for all p G (1, oo) 

with 

p 



pd if p & [2, 00). 

Proof. It is a direct consequence of Sobolev embedding and interpolation between 
LP spaces. □ 

Now we extend the regularity of Q^'^ to general non-smooth kernels. 

Lemma 5.4. Let B be a collision kernel satisfying li2.S\) - (K4\ )- Then, for all p > 1, 
there exist constants C, k and q < p (q only depending on p and d), such that for 
all (5 > 0, and for all measurable function f 

||Q^'+(/,/)L. < c^r'^ll/IUill/IU. + 5||/||li||/IU.. 

Proof. We use a decomposition approach and split the operator Q^'^ as the sum of 
a smooth part and a non-smooth part 



ns ' 



where Qf''^ is given by (15. ip . Then, applying Corollary l5.3l we have for all p G (1, 00), 
there exist q < p, namely (the role of q and p are exchanged here with respect to 
Corollary 15.30 

q=) d+{d-l)p 

1^ ^ if p G [2c/, 00). 

and CregijI^Bf), depending on the the regularization parameter 77 and blowing-up 
polynomially when 77 — > such that 

(5.5) ||gf'+(/,/)L, < ae,(r/,i3f) ii/iui ii/iu,. 
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Now, we need to estimate the remainder Q^^^ = Q^'^ — Qf'"*"- To this aim, we spht 
it as 

Q^s\f, f) = Qi'^if, f) + /) + Q^'^if, f) + •'"(/, /), 

with 



Qf •+(/,/) 
Qf'^(/,/) 



B{\y\,\z\)6{yz)x''{\z\) [x\\y\) - X^ilvD] f f^dydz, 

B{\yl\z\)6{yz)x^{\y\) [x'^d^D - f fidy dz, 

\y\ 



B{\yl\z\)6{yz)x^{\z\)x:;{\y\)Cv 



1 - 



f K dy dz 



B{\yl \z\) 5{y ■ z) x^{\z\) x^i\y\) % (1 - Q f /I dydz. 



On the one hand, we give a first estimate in L^{T>l) applying directly the estimate 
in Lemma [4.11 

(5.6) ligf'+IUi < Ci(i?,S)||/|Ui||/|Ui, a €{1,2,3,4}. 

On the other hand, we treat for instance the operator Qi'^{f,f) and have for a 
fixed V G Vl 



< 



< 



< 



B{\yl\z\)5{yz)xWz\) \xWy\) - X^i\y\)\ \f\\f:\dydz 
B{\v -v'X\v- v'\) 5{{v - ■ {v - v')) 



X \x\\v-v:\)-x^{\v-v:\)\ \f:\dv:)dv' 



LI ||;||l- sup 



BWz\M) dz 



R( 
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< C^{R,B)ri\\j \\j \\L°°- 
Using similar techniques, we prove that for a G {1, 2, 3} 



(5.7) 



\\Q 



R,+ \ 



< CUR.B)v\\fh 



For the fourth term, we have using the cancellation of the collision kernel B at small 



7- 



relative velocities as |f — 
(5.8) WQf^WL^ < Coo(i?,S)r/^||/|Ui||/|Uo 
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Finally, by the Riesz-Thorin interpolation Theorem, from fl5.7ti5.8l) and fl5.6l) . we 
deduce that for p G [1, +00] there exist Cp{R, B) > and jS G (0, 1] such that 

(5.9) IIQf'+(/,/)IU. < CMB)v'\\f\\L^\\f\\Lr. 

To sum up we have obtained for all p > 1 and rj > 0, there exist C > 0, q < p, 
Kq > and P G (0, 1) such that 

WQ'^'HfjnL. < cr-«ii/iu.|i/iu, + r/^ii/iuiii/iu.. 

The conclusion follows by choosing 77 small enough. □ 

Proof of Proposition \5.1[ Now the proof of the propagation of regularity bounds 
is done exactly as in Section 4 and Subsections 5.1 & 5.2] (except for the 
simplification that there is no moments estimates to take care of). □ 

5.2. Entropy — entropy production inequalities. The periodized equation fl3.2ip 
preserves non-negativity, and for a non-negative distribution / one can formally 
compute an H theorem (see [13]): 

jH{f{t)) = -D{f{t)) < 

with 

H{f)= [ f log fdv 

and 

DU) = - I Q^'ifJ) log fdv = ] [ U'f-fh) log B{y,z)dvdydz. 

Then we can state the result which relates the entropy functional H and the 
entropy production functional D: 

Proposition 5.5. We consider the periodized Boltzmann collision operator for some 
truncation parameter R > V^L, and we assume that the collision kernel satisfies 
B > bo \ v — v^l"' , '-f > , for \v — v^,\ < R. Then for any ri,a > there is k E N and 
M, K > (depending only on rj, a, Bq, JjR) such that 

D{f)>KH{f\m^Y+\ m^ = -^ p=[ f dv, 
for any a < f E V-{T>l) with H^^^{T>l) norm hounded by M. 
Remarks: 

1. Note that this is a functional inequality independent of the flow of the Boltz- 
mann equation itself. 
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2. In the case of the classical Boltzmann equation with f G M"^ the entropy 
production functional controls the relative entropy according to the Maxwellian 
equilibrium. Here the equilibrium is a constant, defined by the mass of / divided 
by the total volume of the torus. 

We shall adapt the method developed in [311 Proof of Theorem 2.1]. In the first 
step, we treat the case of a collision kernel B which is uniformly bounded from 
below. In this case we prove the equivalent of the so-called Cercignani conjecture in 
the context of the Boltzmann operator periodized in the velocity space. 

Remark: Note that the assumption R > ^/2L allows to replace in the bound from 
below, when needed, the truncation by the integration over the whole torus. 

Lemma 5.6. Let us consider a collision kernel B which satisfies B > Bq > 0, a 
truncation defined by or liy.16]) . together with R > ^J2L. Then there is 

an explicit constant K such that for any < / G L}{T>l) we have 



Proof. We proceed in several steps. 

Step 1. Since the entropy production functional is monotonous in terms of the 
collision kernel it is no restriction to replace ;B by 1 in the sequel for the estimate 
from below. Moreover it is always possible to bound from below the truncation 
|f — f'l < R and |f — < R (in case we performed the truncation for the fast 
spectral method) by the classical truncation |f — < R. 

Step 2. Using Jensen's inequality on the sphere integration (coming back to the 
classical truncation by the previous remark) and the joint convexity of the function 
(X, Y)^ {X - Y) (logX - logF) on IR+ x M+ we compute 



Let us study more precisely the function G. As it was already observed by Boltzmann 
himself, the function G only depends on f + f ^ and (|f p + |f^p)/2. Moreover, here 
it is also periodic on the torus Vl since / is periodic. It implies (when / is smooth, 
but we can always use mollifications to relax this assumption here) that it in fact 
only depends on t> + t>^. 

Step 3. Let us denote by St the semigroup of the heat equation on L^{Vl) (and for 
brevity we keep the same notation for its semigroup in L^{T>\)). Then the semigroup 



DU)>KHU\m^). 




where F = f f^^ and 
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is compatible with the symmetries in the sense that: 

Stiff.) = (Stf) (Stf.) 

and StG only depends on v + v^, (this follows from a straightforward computation 
using the explicit formula for the Green kernel of St). 

Step 4. Then we have the following computation as in |31j : 

2 



d 
It 



t=0 



St{{F- G) log - {StF - StG) log 1^ 





VF 


VG 




F 


G 



(F + G), 



where V denotes the gradient with respect to € M?'^. 

Then we bound from below the truncation l_B(t,,R)(f*) by the integration over the 
whole torus for (since R is large enough), and we compute 



d 
'di 



D{Stf) 



t=o 





VF 


VG 


/ 


F 


G 


-I M 


\f-( 




VF 


VG 


j 


F 


G 



(F + G) dv dv^ 



F 
G 



(F + G) dvdv^. 



We deduce by the semigroup property that for alH > 



WSfF VStG 



StF 



StG 



and therefore 



D{f)> 



+ 00 



VStF VStG 



StF StG 
Step 5. We now use the fact that the operator 

P : 



{StF + StG) dvdv. 



{StF + StG)dvdv^ dt. 



{A,B) 

is bounded from x to V^. Hence 

VStF VStG 



{A-B) 



StF 



StG 



2 


PVStF 




StF 
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since PVStG = V^G — V^^G = from the fact that G only depends on v + v^,. We 
deduce 

{StfSth + StG)dvdv, 1 dt 



D{f) >C / 

Jo \JVlxVl 

and thus (dropping the term SfG) 



Stf V Stf 



D{f) >G [ 

Jo \JVlxVl 



Stf V Stf 



Stf Stf^ dv dVi, dt. 



Step 6. From now on the proof departs shghtly more from ^31j: it is simpler since we 
are in the torus and we have more symmetries. Let us show the following functional 
inequality: for any smooth non-negative function h, 



L 



h 



h 



where 



nh\g) 



h hi, dv dVi, > G2 I{h\mc 



dv. 









Vv log - 
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The proof only amounts to Jensen's inequality on the variable f^: since J^, hdv = p, 









1 

IVlxVl 


h 





h hi, dv dVi, > 



C 



Then as 



we deduce 



P Jv^ 
Wh^dv^ = 



6Cl 



Vyh f'Vv^h 



h 



h 



h dv. 



Iv, 



h 



h 



Step 7. So far we have proved 



hK> G 



V„h 



hdv = C I{h\moo). 



r+00 

D{f) > D{f) > G, / I{Stf\ 
Jo 



rrioo) dt. 



Then a trivial computation shows that 

d 



= -I{Stf\mc 

Moreover from the explicit formula for St we have 
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and thus we finally obtain 

D{f) > D{f) > Cs iH{Sof\m^) - 0) > CsH{f\m^). 

□ 

Now we are ready to prove Proposition [531 Since we deal with a bounded velocity 
domain we do not care about possible decay of the collision kernel at large relative 
velocity (as for soft potentials) and the only cancellation we have to treat is for zero 
relative velocities. 

Proof of Proposition 15.51 We only mention the difference as compared to the previ- 
ous proof. 

The reduction to a collision kernel uniformly bounded from below studied in 
Lemma [5.61 is done as in ^31j [Theorem 4.1]: one write for some small 6 > 

B{\v-v4) > 6^ {Bo-lBioM\^-^A)) 
where Bq > bo > and we deduce 

D{f) > 5^ {DoU) - DsU)) 
where Dq is the entropy production functional corresponding to -Bq, and 

Then we have the following Lemma, which is proved exactly similarly as [3T1 
Theorem 4.2]. It is even simpler since Maxwellians are replaced by constant functions 
and the study of the tail is not needed (we omit the proof for brevity). 

Lemma 5.7. For any e G (0, 1) and a > 0, there are constants of smoothness k, M 
and some corresponding constant Cdiag > such that 

D5<C,,agHif\m^Y-'S''/^ 

for any a < f E L^(Vl) with Hp^^iVi) norm hounded by M. 

But since we have 

DoU)>CHU\m^) 

from Lemma 15. 6[ it is straightforward to get the result by choosing correctly the 
parameter 5. 

□ 

Now we can proceed to the proof of Theorem 13. 1[ 
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5.3. Proof of the global in time stability. In this subsection we shall turn to 
the question of obtaining uniform bounds as well as global existence, in order to 
conclude the proof of Theorem 13.11 Indeed in Section 4 the smallness assumption 
on the truncation parameter e a priori depends on T and could go to as T goes 
to infinity, since it depends on regularity bounds growing exponentially in times. 
In order to overcome this difficulty, we shall combine the following arguments: 

• for the unperturbed problem (13.211) we have a Liapunov structure and the 
solution converges to a unique prescribed equilibrium from the regularity 
and entropy production studies; 

• the equilibrium distribution of the unperturbed problem fl3.2ip (that is the 
constant functions on the torus) are also equilibrium distribution of the per- 
turbed problem (13.171) : 

• by taking the size of the perturbation small enough (measured in terms of 
e) it is possible to construct a solution to the perturbed problem on an arbi- 
trarily large time interval [0,T], on which moreover the perturbed solution 
remains close to the unperturbed solution (I3.2ip . say in Sobolev norms; 

• finally the constant equilibrium functions are non-linearly stable for the per- 
turbed problem, with a stability domain independent on the size of the per- 
turbation. 

Hence we shall deduce that as soon as the time for which the perturbed solution 
departs from the unperturbed solution is larger that the time-scale of relaxation 
to equilibrium for the unperturbed problem (I3.2ip . the perturbed solution shall be 
trapped by the stability domain of the equilibrium before instability due to the 
perturbation can develop. Let us formalize these arguments in a proof: 

Proof of the global stability and asymptotic behavior in Theorem \3.1[ Let us consider 
some initial datum < /o G which belongs to Hp^^iVi), and some smooth 
balanced perturbations /o,e of it. These perturbations have the same mass, and 
therefore are corresponding to the same equilibrium (this is the reason for this as- 
sumption). 

Since the smooth balanced perturbation preserves the constant equilibrium of 
(I3.2ip and is nonlinearly locally stable in any Hp^j.(Vi), we fix a 77 > such that 
the constant distribution moo = p/I'^lI associated to the mass p of /o on the torus 
has attraction domain with size rj in H'^^^{T>l) for the perturbed problem (I3.17p . 

On the one hand, we show that there exists a unique solution f{t) to (13.211) and 
from Propositon 15.11 we obtain uniform regularity bounds for alH > 



ll/Wlk,^.<c. 
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Moreover, from Proposition 15.51 there exists a time Tq such that the solution f{t) 
is ri/2-close to the equihbrium in H^^^iVi) (p < k) for t >Tq (using the Csiszar- 
KuUback inequahty in the torus, see Theorem 1] for instance): 

\\f{t)-m^\\ni^^, < V2. 

On the other hand, applying Propostion 14.61 with T = Tq, we prove that there 
exists e, which only depends on the Hp^^iV^) and L^iV^) norms of /o such that for 
all e such that < e < i, there exists a unique smooth solution f^(t) G Hp^^iVi) to 
(IXTTj) on [0, To] such that 

VtG [0,To], mmm,,. < C(To), 

and (for any p < k) 

\\f{t) - feimn^,^,. < VToie), 

where f{t) is solution to (13.211) and <^To(£^) goes to zero when e goes to zero. 

Then, we fix a perturbation parameter i small enough such that for e G (0, e) the 
perturbed solution satisfies 

Finally, at time Tq the perturbed solution belongs to the stability domain of the 
constant distribution for the perturbed problem and it is trapped. 

Therefore there exists a unique global smooth solution /g, which is uniformly 
bounded for all t > and such that for any p < k 

\\m\\H^,.. < max(C(To),C + r^). 

This achieves the proof of (i), (ii), {Hi) and {iv). 

□ 



6. Application: stability and convergence of spectral methods 
In this section we consider the following spectral approximation of (13.211) 

= VNQ^ifN, /at), 

where Vn denotes the orthogonal projection on Pjv in L'^{Vl) (the space of trigono- 
metric polynomials with degree less at most N in each direction). 
The goal of this section is to prove the following theorem: 

Theorem 6.1. Consider any non negative initial datum fo G Hp^.^{T>L), with k > 
d/2, which is not zero everywhere. Then there exists Nq E N (depending on the 
mass and Hp^^{T>L) norm of f) such that for all N > Nq: 
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(i) there is a unique global solution fj^ = /Ar(t, ■) to the following problem 

'Pn Q^ifN, In), 



dfN _ ^ 



(6.1) 



dt 

fN{t = 0) = VNfo; 



{a) for any p < k, there exists C > such that 

Vt>0, \\fNit,-)\\H^^^ < C- 



(Hi) this solution is everywhere positive for time large enough, and the mass of 
its negative values can be made uniformly (in times) L°° small as N ^ oo; 

(iv) this solution /jv converges to f{t) the solution to ^3.21\) with the spectral 

accuracy, uniformly in time; 
{v) this solution converges exponentially fast to a constant solution on the torus 
prescribed by the mass conservation law. 



To prove Theorem 16. we want to apply Theorem 13.11 with the perturbation 

PnUn) '■= Vn Q^ifN, fN) — Q^ifN, fw), 

which preserves the mass: 

/ P^{f^)dv=[ (rNQ^'ifN^fN) - Q^'ifNjN)) dv = 0, 

In the next Lemma, we prove a consistency and smoothness result for this ap- 
proximation. 

Lemma 6.2. Consider a non negative function f G H^^^{T>i), with k > d/2, which 
is not zero everywhere. Then, there exists C > depending only on the collision 
kernel B and the truncation such that for all p G [0, k] we have 

(6.2) WPSmnl.. < C||/|Ui||/||h,^... 

Moreover, for all p E [0,k] 



(6.3) \\P^{f)\U.. < C\\fU."'""^- 



]\fk-p 



Proof. First, we split the operator P/v as 

WP^iflWHU < \\Q''ifJ)\\H^,.r + \\PNQ\f,f)\\H^,^., 

As in the proof of Lemma [4.21 we get that for all p G [0, A;] 
(6.4) IIQ^(/,/)llk. < C,{R,B) Il/IILII/Ilk, 
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Concerning the interpolation error estimate, the following result holds. If m G 
-^per(^i) f*^^ some p > I, then 

C 

(6-5) \\u - VnuWh^^^^ < Mm,^^- 

Then, taking p = A; in the latter inequality and from fl6.4l) we obtain 
(6.6) WV^Q^UWh.^^ < \\Q\fjm^^^^ < C,{R,B) WfWl 11/11^. 
Gathering (16.41) and (16. 6p . we finally get 

< C,{R,B) ll/IUi 11/11^. 

Moreover, using again the error estimate (16.51) . it yields 

\mfj)\\HL. 



per 



^per 



]S[k-p 



□ 



Let us now perform a linearized study of the perturbed equation (16.11) by classical 
Fourier-basis decomposition. The only equilibrium distributions of the equation 
(I3.2ip are the constant, prescribed by the mass conservation. Let us consider the 
linearized version of the perturbed equation (16. ip around such a constant moo'- 

where 

/:^.«(/) = P^ [Q«(/,i) + g«(i,/)]. 

Let us prove the following lemma: 

Lemma 6.3. The operator is bounded and self-adjoint in L?'{T>l). Moreover 
it is non-negative, its null space is given by the constant functions, and it has a 
spectral gap A > 0. As a consequence, the constant are nonlinearly locally stable in 
any Hp^^{T>L) for the equation (16. ip . with a stability domain independent on N . 

Proof. The boundedness is trivial. Then, the periodized operator is translation 
invariant, which implies that the Fourier modes 
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for k G Z*^ are trivially eigenfunctions of This provides a complete orthonormal 

eigenbasis in Li^iVi). A trivial computation yields 

>C^''^(ejk) = Ofe ejfc l|fe|<jv with Ofc := - / [l + ek{y + z)-ek{y)-ek{z)\B{y,z)dydz. 

Jcr 

In particular we deduce that C^'^ = C^'^Vn, and the self-adjointness comes from 
the following identity obtained by the usual changes of variables: 

X [Pnq' + Vnq'^ - Vn9 - VNQi] B{y, z) dy dz dv. 
Another formula for is readily deduced from the previous representation: 

Ofe = -7 / {{ek)' + {ekYi,- (ek) - {ek)*\'^ I3{y,z)dydzdv. 

One sees from the second representation that = a-k < for any k e Z*^, and 
from the first representation and Lebesgue theorem it is easily seen that for — > 00 
the coefficients ak converge as 

— Ooo = - / B{y,z)dydz e (-00, 0). 
I'^Hoo Jcr 

Hence we deduce that a^. G [aoo,0] for any k G Z'^, with asymptotic convergence 
towards Ooo for large k. Moreover the null space can be computed: for some smooth 
periodic function 0, the equation 

0' + 01 - - 0, = 

implies that the third-order derivative of is zero, and the periodicity then imposes 
that it is a constant. Thus the null space is spanned by Cq. It concludes the proof 
of the existence of a spectral gap 

Xn :=min{|afc|, ke [|-A^,A^|]^ k O} > 

which is uniformly bounded from below as — > +00, since 

Ajv — ^ Aqo '■— min {\ak\, keZ"^, /c 7^ 0} > 0. 

The non-linear stability in comes from the fact that for the perturbation h — 
f — TTioo, we have the following control on the bilinear part: 

\\VNQ''ih,h)\\H^^^<CB,R\\h\\l2 

for some given constant Cq^r > independent on N. 

Finally using the eigenbasis of the Fourier modes, a similar study can be performed 
in any Sobolev space i7pg^('DL). □ 



ak 



ANALYSIS OF SPECTRAL METHODS FOR THE BOLTZMANN EQUATION 



39 



Remark: Exact computations could be made for particular physical collision ker- 
nels B (in a similar way as the computation of the kernel modes (3{l,m) in [TTl 

mm)- 

6.1. Proof of Theorem 16.11 Consider the numerical solution given by solving 
(16. ip . We can formulate the problem as a perturbation of the truncated Boltzmann 
equation. Indeed setting 

PnUn) = -(/rf - Vn)Q''Un, In), 

the problem (16.11) can be written as 

dfN 



dt 



Then, applying Lemma 16. 2[ the perturbation P/v satisfies the assumptions of Theo- 
rem [SHJ Moreover since /o G H^^^iVL), we have straightforwardly 

il/A^(0) linger ^ ll/olln^er' ll/jv(0) " /ollff^er ^ 0" 

Hence, we can directly apply Theorem 13.11 to the perturbation P/v, which proves 
that there exists A^o large enough and only depending on /q, the kernel B and 
the truncation, such that for all N > Nq, the perturbed system admits a unique 
uniformly smooth solution, which converges to a constant, and satisfies all the points 
in Theorem 13. 1[ 
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